clc,clear
U=230;
R=100;
L=265e-3;
w=2*50*pi;
Z=(R^2+(w*L)^2)^0.5;
fi=atan(w*L/R);
A=2*pi/3; %ugao paljenja u radijanima
tau=L/R;

b=270*pi/180;
uslov=1;
while uslov
    f=2^0.5*U/Z*(sin(b-fi)-exp((A-b)/(w*tau))*sin(A-fi));
    f_p=2^0.5*U/Z*(cos(b-fi)+1/(w*tau)*sin(A-fi)*exp((A-b)/(w*tau)));
    
    b=b-f/f_p;
    if abs(f)<1e-4
        uslov=0;
    end
end

vremeVodjenja=b-A %vreme vodjenja tiristora u radijanima
b=b*180/pi %ugao gasenja tiristora u stepenima